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Present expansion stage of the universe is believed to be mainly governed by the cosmological con- 
stant, coUisionless dark matter and baryonic matter. The latter two components are often modeled 
as zero-pressure fluids. In our previous work we have shown that to the second-order cosmological 
perturbations, the relativistic equations of the zero-pressure, irrotational, multi-component fluids 
in a spatially near flat background effectively coincide with the Newtonian equations. As the New- 
tonian equations only have quadratic order nonlinearity, it is practically interesting to derive the 
potential third-order perturbation terms in general relativistic treatment which correspond to pure 
general relativistic corrections. In our previous work we have shown that even in a single compo- 
nent fluid there exists substantial amount of pure relativistic third-order correction terms. We have, 
however, shown that those correction terms are independent of the horizon scale, and are quite 
small (~ 5 X 10~^ smaller compared with the relativistic/Newtonian second-order terms) due to the 
weak level anisotropy of the cosmic microwave background radiation. Here, we present pure gen- 
eral relativistic correction terms appearing in the third-order perturbations of the multi-component 
zero-pressure fluids. As a result we show that, as in a single component situation, the third-order 
correction terms are quite small (~ 5 x 10~^ smaller) in the context of the evolution of our universe. 
Still, there do exist pure general relativistic correction terms in third-order perturbations which 
could potentially become important in future development of precision cosmology. We include the 
cosmological constant in all our analyses. 

PACS numbers: PACS numbers: 04.50.+h, 04.62. +v, 98.80.-k 
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References [l^ 

I. INTRODUCTION 

Recently, we have been presenting a series of work based on our theoretical study of relativistic nonlinear cosmo- 
logical perturbations SB- We have shown that to the second-order perturbations, general relativistic equations 
of a zero-pressure, irrotational fluid in a spatially flat background have exact Newtonian correspondence except for 
the presence of the gravitational wave contributions P, In an accompanying paper [5] we have relaxed all the 
assumptions we made in the second-order perturbations in fT, '3| , and have derived pure general relativistic effects 
from the pressures, rotation, spatial curvature, and multi-component. In that work, we have shown that except for the 
multi-component situation, relaxing any of the other three assumptions leads to pure general relativistic correction 
terms appearing in the second order. Pressures are intrinsically general relativistic even in the background and the 
linear-order perturbations. The presence of background curvature leads to first non- vanishing relativistic correction 
terms appearing in the second order. The rotational perturbations generally lead to relativistic correction terms which 
become negligible in the small-scale (sub- horizon scale) limit, thus having relativistic/Newtonian correspondence in 
that limit. In the case of zero-pressure, irrotational multi-component fluids in a flat background, effectively we have 
exact relativistic/Newtonian correspondence even in the multi-component situation; this will be summarized in the 
later section of this work. 

The relativistic/Newtonian correspondence in the background world model was known in the zero- pressure medium 
by the work of Fricdmann in 1922 [g^ in the context of Einstein's gravity, and by the work of Milne and McCrea in 1934 
[7| in the context of Newton's gravity; the latter Newtonian derivation is later known to be a guided one by the already 
derived Einstein gravity result Q. In the case of linear perturbations, the relativistic/Newtonian correspondence was 
also known in the zero-pressure medium by the work of Lifshitz in 1946 Q in the context of Einstein's gravity, and 
by the work of Bonnor in 1957 in the context of Newton's gravity. The fully nonlinear perturbation equations 
in the context of Newtonian cosmology in a zero-pressure medium were known in a textbook by Peebles in 1980 
[llj . The Einstein's gravity counterpart of the nonlinearly perturbed cosmological medium, especially the continued 
relativistic/Newtonian correspondence even to the second order, was first shown only recently in our works in [1, 2]. 

Meanwhile, in we derived pure general relativistic correction terms appearing in the third-order single compo- 
nent, zero-pressure, irrotational fluid in a flat background. Thus, now it is a natural step to find out the potential 
third-order pure general relativistic correction terms appearing in the multi-component, zero-pressure, irrotational 
fluids in a flat background. As the Newtonian system has only quadratic nonlinearity even in the multi-component 
situation, see Sec. II of Q, any nonvanishing third-order terms can be regarded as pure general relativistic corrections. 
The situation is also practically important because current stage of the universe is supposed to be dominated by two 
zero-pressure components and the cosmological constant. We will include the effect of cosmological constant in all 
our analyses and equations in this work which is also true in our previous works in [l], [2, H, 0) IE] ■ 

In Sec. |TT] we present the metric and fluid quantities perturbed to the third order which will be required in our 
calculation. We present fluid quantities for most general fluids with pressures and stresses which will turn out to be 
important even in the zero-pressure situation in our main analysis. As in the single component case in under our 
proper choice of variables and gauges we do not need third-order perturbations of the connection or curvature tensor. 
In Sec. mil we summarize the relativistic equations to the second order, and their effective correspondence with the 
Newtonian ones even in the multi-component situation. In Sec. llVl we derive the general relativistic third-order terms 
and present equations in the context of Newtonian gravity with pure general relativistic corrections. We compare 
our equations in the multi-component case with the previously derived ones in a single component. Section fVl is a 
discussion. We often set c = 1, but recover c in the Newtonian context presentation. 

II. THIRD-ORDER PERTURBATIONS 

A. The covariant and ADM equations 

In the following we summarize the basic sets of covariant equations and ADM equations we need in our analysis. 
These equations, except for the covariant equations of individual component, are presented Sec. II of [l|l; notations can 
be found in that work. For original studies of the covariant and the ADM equations, see [l^l, and [lj|, respectively. 
Although, we will use the ADM equations in our calculation, the covariant equations show another aspects of the 
same fully nonlinear system of Einstein's equation. 
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The energy-momentum tenor of a fluid can be decomposed into fluid quantities as 

fab = flUaUb + P {UaUb + Qab) + QaUb + qbUa + TTab- (1) 

Without losing generality, we take the energy frame, thus set Qa = 0. This decomposition is valid even in the multiple 

component fluids; in such a case the above fluid quantities can be regarded as collective fluid quantities. In the 
multi-component case we introduce the energy-momentum tensor and fluid quantities of individual component as 

fab = '^T{j)ab, 
j 

T{i)ab = ij-{i)U{i)aU(i)b + P{i) {u{i)aU{i)b + 9ab) + T^{i)ab, (2) 

where, without losing generality, we also took the energy-frame condition for each component, thus set ^(j)^ = 0. For 
interactions among components we introduce 

In a single component situation, taking the energy-frame (thus, setting Qa = 0), the energy conservation equation, 
the momentum conservation equation, and the Raychaudhury equation are 

A+(M+P)^+^"'''^a6 = 0, (4) 

(/i +p)aa + hi {p,b + ^6%) = 0, (5) 

^ + 1^2 _ -a^^ ^ -ab-^^ _ -ab-^^ ^ + 3^) _ A = Q. (6) 

In the multi-component case, taking the energy-frame for individual component (thus, setting ^(j)^ = 0) we have 

+ +P{l)) ^(i) + T^(i)^{i)ab = -u'li)I{i)a, (7) 

(Aw + P{i)) a(i)a + (Pii),b + ^(i)b-c) = ~h{i)''aki)b' (8) 

^(^) + l^l) - «W;a + ^ii)^(^)ab " ^i-)^ir)ab " 47rG [fl - 3p - 2fabUy)u\,)) - A = 0. (9) 



In order to handle cosmological perturbations to the third order, we also need the momentum constraint equation for 
a collective fluid 

hab - Sr\ + ''^ + {Qjab + ^ab) = 0. (10) 

In the ADM formulation, the energy conservation, momentum conservation, and trace of ADM propagation equa- 
tions are 



(11) 



J^,oN-^ - J^-.pN^N-^ - JpN^.,^N-^ - KJ^ + EN^^N-^ + S^^ + S^N^pN-^ = 0, (12) 

K^qN-^ - K^^N°'N-^ + N'-'^^N-^ - K°'^K„i3 - \k'^ - 4ttG {E + S)+A = 0. (13) 

In the multi-component case we have 

E^i^^oN-' - E^,)^^N'^N-^ - K {e(^ + ^-S(^ - S'^^K^p + N'^ (^'-^w).^ = (^^(^)o " , (14) 

The momentum constraint equation is 

^f./3-^^,« = 87rGJ„. (16) 
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The above sets of equations are only parts of the covariant and the ADM equations; for complete sets, see [U, [H, 
We will show that, by taking proper choice of gauges, the scalar-type perturbations to the third order can 
be derived from either of the above sets of equations. We will present the derivation based on the ADM equations, 
because the covariant formalism often requires lengthier calculation in our particular case; of course, the covariant 
equations also give the same result. As we will show, however, we use both formulations simultaneously depending 
on the convenience. 

B. Metric 

Our metric convention is the same as in [l| 

300 = -a^ (1 + 1A) , goa = -a^S„, ~g^p = (5^ + SCa/j) , (17) 

(3) 

where tensor indices of and Ca.fi are based on (7^^. To the third order, the inverse metric becomes 

5°° = -4 (1 - 2^ + - - + iAB^Ba + 2B°'B'^Cai3) , 

= --^ [B" - - 2B'^C;^ + (4^2 - Bf^Bp) + {AB^^ + B'*C% , 

~g°^i^ = ^ (^5(3)"/3 _ 2C"/3 _ 5"s/3 ^c^c^-i + 2AB°'B^ + 2B°'B''&^ + 2B>^B''C^ - 8C;^Cf C*) . (18) 

In order to derive perturbation equations to the third order, we need the connection only to the second-order. These 
are presented in Eq. (52) of 

The ADM metric variables follow from Eq. (2) of [l| as 

N^a(^l + A-^A^ + ^B^^Ba + ^A' - ^AB^'Ba - B'^B^C^^ , 
7V„ = -a^Ba, ^ -B" + 2B^C'^ - AB^^C^Cj, 

Kp = a' (.gS + 2^^"/3) , h'^^^^ (5^'^"'' - 2C"'^ + ^C?,Cf^ - SC^Cf C^^) , (19) 
where tensor index of is based on h^fi as the metric; h°^^ is an inverse metric of h^fi ■ 

C. Fluid quantities 

We introduce perturbations of fluid quantities as 

/i = ^ + (5^, p^p+Sp, Ua=aVa, TTap = (20) 
(3) 

where tensor indices of Wq and IIq^ are based on g^p- The above fluid quantities can be regarded as collective fluid 
quantities in the case of multi-component fluids, see below Eq. (|26p . Although wc will consider zero-pressure fluids, 
it is important to keep the perturbed pressure (Sp) and anisotropic stress (Hq,^) because the collective pressure and 
anisotropic stress do not vanish to nonlinear order in the multi-component fluids even in the zero-pressure case, see 
Eq. (j44p . In any case, for later convenience, in this section we will present fluid quantities for most general fluids. 
Components of the four-vector Ua are 



Ua = aVa 

uq = —a 
a 



1 + A - i + i («« + B") [va + Ba) + \a^ + \a {v^v^ ~ B^B^) - C^p (^'" + B") {v^ - 
+ - AS" - 2C"'3 {vp + Bp) + ^A^B" + 2AB^^C$ + ^B" {v^vp - B^Bp) + AC$C^ {v'' + B^) 



2 

■ (21) 



l-A+^A^ + \ {v'^va - B'-Ba) - \a^ ~ \a {v'-va - SB^^B^) - C^p {v'^v^ - B'^B^) 
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In instead of Va, we used Va defined as 



= V". (22) 



Thus, we have 



= V^a - S„ + ABo, + 2Co,pVl^ - \a^B^ ~ B^Vf3 {^V" - B^^ . (23) 

Components of iTab are 

= a^n,^, noa = -a'n„^ [v'^ + S'' + Aw'' - 2C^^ (i;^ + B^)] , noa = a'n„^ (w" + B") {v^ + B^) . (24) 
From TT^ = we have 

jja _ 2C7"''n„;3 + (4C^C'3^ - v'^v^) n^p = 0. (25) 
To the thnd order, the energy- momentum tensor becomes 

f = -/i ~Sfi~{fi + p) v" («„ + B^) + ifi+p) [AB^ + 2C„/5 {v" + B^)] v" 

- {Sfi + Sp) (w„ + S„) - n,;,,;" (v'' + B^) , 

- {Sfi + Sp) Ava - (Avl^ + 2C'^^v^) , 

2™ = + ,5pJ| + n;^ + (^ + p) („" + 5") _ 2C'"'^n^^ - (/i + p) [AB" + 2C"^ {v^ + B^)] vp 

+ {5^l + Sp) (u" +B'')vi3+ Up^ (S"w^ + 4Q"C^*') . (26) 

Above relations of four-vectors, energy-momentum tensor, and fluid quantities are derived for a single component 
fluid. However, these are also valid as the collective component in the case of multi-component fluids. We can 
easily see that, by replacing these quantities with the ones of individual component, the same relations are valid for 
individual component as well. That is, for an individual component, say i-th component, by replacing 

Tab, A; P: Ua, TTab-, P> ^^i, Sp, Va, Va, ^afi, (27) 

with 

T{i)ab, P(i), U(i)a, T^{i)ab, : 7 5P(i), 11(^)0^, (28) 

Eqs. are valid for the «-th component. 

From Eqs. pjl. ipS)) we can express the collective fluid quantities in terms of the individual one. To the background 
order we have 

3 3 



To the third order in perturbations, we can show 

S^x + [in +p)v'' + {5,1 + 5p) + W^pv^] Va-2{fi + p) C^^ffV^v^ 



2(/^0)+PO))Ca/3^^0)^^f,)}, 



3 

h+\ [(m + p)v'' + {5fi + 5p) v'' + Yl'^v'^] v^-^iii + p) Ca0v"v'^ 

^ C \ r 1 2 ^ 



{lJ, + p)Va + 



Sn + Sp+ -{lJ,+p) V^Vff 



■^/3 - 2 (/"+P)^a^''^/3 



+ 



2 



^uh - 2 im + P(3)) '^^s^uA3) } f • (30) 



Thus, we have 

Sf^ = Y, {'^^o) + (^0) +^'0)) ^0) (%)« - ^«) 

<5p = ^ {SPU) + I im +PU)) ^U) (^U)a - Va) 
3 

+1 [i^m + ^P(3)) ^{3) + - 2 +p{3)) c^v^j)] (?^(,)a - v^) }, 

3 I 

j I L J 

r 1 2 

+ [i^m + ^pu)) ^(3) - ^u"y\ - - 3 (^0) +Pi3)) (^0)7 - ^7) 



A" 



(31) 



The ADM fluid quantities are introduced as 



E = haUbT"''', Ja = —ntTj^: Sap = Tap, S = h^^Sap, Sal3 = Saff — -^hapS, 



(32) 



where tensor indices of Ja and Sap are based on hap as the metric. The four-vector ha is a normal frame four-vector 
with ha = 0- Thus by setting z;^ = we have Ua = ha, and Eq. (|2ip gives 



ha = 0, 



no 



a[ l + A--A'^ + -B"B. 



Ia^ - lAB'^Ba - CapB°'B'^ 



re" = i (^B" - AB°' - 2C"^B/3 + ^A^B"' + 2ABl^Cj^ - ^B^Bf^Bp + AC'j^C^B'^ 
l(^l-A+^A^- \B'^Ba - ^A" + ^AB<-Ba + CapB^^B^^ . 



Using Eqs. dS!]),®, Eq. ([22]) gives 



E = fi + Sj.l+ (/i + p) V^Va + [5^ + 5p) V^'Va - 2 {fl + p) Capv"v'^ + Uapv'^V 



Ja — a 



{fl + p)Va + {SfJ, + Sp) Va + UapV^ + + P) VaV'^Vp - 2IlapC^v'< 



Sap = pgl^^ + SpQap + n^/J + 'ipCaP + {fJ. + p) VaVp + 2SpCaP + {S^i + 6p) VaVp 
S^ip+i6p+ (/i + p) V^Va + {Sfl + 5p) v'^Va + ^ [HaP - 2 {fl + p) Cap] , 

Sap ^0.^1 IIq/J + (M +P) (vaVp - o^i^/S^'^^T ) + (^M + Sp) VaVp - ^ i^J. + p) Capv'^Vj 



IdS [(^M + Sp) v-^v^ + n^sv^v' -2{fi+p) C^sv-'v'] 



(33) 



(34) 



The individual ADM fluid quantities can be found by replacing 

E, Ja, SaP, S, Sap, (35) 

with 

J(i)a, S(i)ap, S(i), S(i)ap, (36) 

and similarly for the fluid quantities and the energy- momentum tensor as in Eqs. (P7)) . (P5|) . From Eq. ^ we have 

£■ = ^ E(j), Ja^^ Ji^j)a, Sap = ^ S(^j)ap, -5" = ^ S'(j) , Sap = ^ S(j)ap- (37) 



D. Decomposition 

We decompose the metric into three perturbation types [l^ 

A = a, Ba= P.a + B^:^ , Cap = (^ffg + 7,a|/3 + CgJ^) -|- C ^ *j , 



(38) 



where superscripts (v) and (t) indicate the transverse vector-type, and transverse-tracefree tensor-type perturbations, 
respectively. Only to the linear-order perturbations in the homogeneous-isotropic background, these three types of 
perturbations decouple and evolve independently. We introduce 



which are spatially gauge-invariant combinations to the linear order We set 

K = -3H + K. 



(39) 
(40) 
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By using n we can avoid third-order expansion of the trace of extrinsic curvature K . Identifying k, with Newtonian 
velocity variable later will be an important step in our analysis, see Eqs. (|63p . (ll0ip . 
For the fluid quantities we decompose 

The perturbed fluid velocity variables v and Va'' subtly differ from the ones introduced in Q, see 0. 



E. Zero-pressure irrotational fluids in the comoving gauge 

The zero-pressure condition sets 

Pit) = 0, EE = n(i)„^. 

The comoving gauge condition {v = 0) and irrotational condition {va^ = 0) give 

Va = 0. 

Under these conditions, Eq. (|3ip becomes 

j 



(42) 



(43) 



(44) 



Notice that for the nonlinear perturbations, the collective pressure bj) and anisotropic stress IIq,^ no longer vanish 
even for zero-pressure fluids; also, the collective energy density (5/x is no longer a simple sum of individual component. 
For the ADM fluid quantities, Eq. ([34]) becomes 



E = fi + Sfj., Ja =0, Sap = {Spg^al + ^a/j + ^SpCapj , S = 3(5p, Sap = a^^ap- 
For the individual component, Eq. (j34p gives 

'S'(i)q;3 = {fJ'{i)V{i)aV{i)p + <5Ai(i)W(j)a^'(i)/3) 7 



(45) 



S(i)aP — 0^ 



1 



(^W(i)a'y(«)/3 - ^9apvjif{i)'i j + SfI(^)V^i)aV(,)p - ^f^{t)CapvJ^^V(^^)-, 



(46) 
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III. SECOND-ORDER EQUATIONS: SUMMARY 



In this section we summarize our previous work of second-order perturbations in zero-pressure, irrotational, but 
multi-component fluids @. The results show that, to the second order, effectively the relativistic equations coincide 
with the Newtonian ones even in the multi-component situation. This provides a reason to go to the third order 
in relativistic perturbation in order to find out pure general relativistic deviations from Newton's theory. Since the 
zero-pressure Newtonian perturbation equations have only quadratic order nonlinearity, any non- vanishing third-order 
terms in relativistic analysis can be regarded as a pure general relativistic correction terms. In Q we presented such 
third-order correction terms in a single component case. In the next section we will derive the third-order correction 
terms appearing in the multi-component situation. As the third-order analysis closely follows the case of second-order, 
in the following we will derive the second-order equations directly from the ADM equations presented in Sec. Ill Al 

We consider zero-pressure multi- fluids, thus 



= 0, = = n(i)„/3. 



(47) 



The irrotational condition and the temporal comoving gauge condition lead to Va = 0, thus we have Jq, — 0; as the 
spatial gauge condition we take 7 = 0, thus /3 ~ x/a. To the second-order perturbations, Eq. IT^ gives 



(48) 



thus 



N = a- 



1 a 

2 



/3 



(49) 



Using Eqs. gSD,®, Eqs. (HI]) , (UnD , mi) , (HH) give 



P 

0)1/3 



k + 2Hk - AttGqS 



3iJ + ^tiGq + (? 



K + C — 

a 



i [a (I + %)) ■ = --^ (ww/3X'^) „ + - ^ (^'Wa^'f,)) 
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c 

"2a 



1/3 



7 /3 

^(i)%)l/3 



(50) 



(51) 



(52) 



(53) 
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From Eqs. (|50l) , (|5T]) , and Eqs. (|5T])-(|53l) we can derive 



1 



a' (S + 4'5,"X'")1 - ^T^GgS (1 + 5) - -4'^,«X'" + 

\ / J Qj O 



1 / A 



\ Q, / J Qj O 



(54) 



1 / A 



(k + (/j)'" + 2 (k - V3) 



(55) 



We have recovered c using 

[G^?] = ^-^ b] = [/i] - [^^c2], H = [i?^')] = i, [v] = i, 

[^] = = [7] = [Cl^l] = K] = = W = = 1, [x] = [«] = 

where i?^'^) = 6K with normahzed ^ = or ±1. 
Without the rotational mode, we introduce 

Since we used the comoving gauge condition, to the linear order, we have 

where = v — x/a and vaw = i'(i) — x/a to the linear order. To the linear order, Newtonian velocity perturbation 
variables are introduced as [l7| 



(56) 



(57) 



(58) 



thus 



Dimensions are 



Thus, we have 



V(j) = -Vt)(i) = - (u, - u) 



[i;] = = 1, [u] = [u,] = = [m,] = [c] = L/T. 



Y: ^A-V. = E f . [(u, - u) V . (u, - u)] 



(59) 
(60) 
(61) 



(62) 



A. Newtonian correspondence 

We assume a flat background, thus = g Thus^ we have = to the linear order even in the presence of the 
cosmological constant and multiple components [l3,[i3|, see Eq. ((99)) . To the second order, we identify the Newtonian 
perturbation variables S, 6i, u, and as 

K„ = -iv-u, u = Vu, v(j)„ = i (ui - u) , S = 6y, 6i = S^i)y. (63) 
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Comparison of the consequent relativistic equations with the Newtonian equations will show apparently why these 
identifications are the proper ones. Examination of Eqs. (|50p - (|55l) shows that we need x{— Xv) only to the linear 
order. From the momentum constraint equation in Eq. (fT6|) . see Eq. (197) of [l|, we have 



a 

thus 



Cz^Xv, (64) 



Xv = -u, (65) 



to the linear order. 

Equations ((50l) - ((55)) become 



j + iv • u = --V • (,5u) + i/^ ^1 |i|u, - up + SA-iy ■ [(u, - u) V • (u, - u)] j , 

-V • (u + Hn) + AtiGqS = -^V • (u • Vu) - C'^*)"'^ + (jiA 

a \a ' ) 

+ E 7^ {^^Ge ~ c^^) |u, - up + {xI^Gq - c^^) A-iy . [(u, - u) V . (u, - u)]| 
+ • u, - _iv • ((5,u,) + - [2(^V • (u,; - u) - (u, - u) • + 2 « - u'^f 



a a a 



+H^\n. _ up + 3i/ ^ ^1 |i|u, - up + A-^V • [(u, - u) V • (u, - u)]| 



-V • (u^ + Hui) + AttGq5 = -4tV ■ (ui • Vu,) - C^*)"'^ ( -u^ 



*l/3 + Cfjj 



+A^GY,q,^ |i|u, - up + SA-iy • [(u, - u) V • (u, - u)]| , 

1 (aH) - A^Gg6 = [aV ■ {5n)] + • (u • Vu) + C7(*)"^ {^^u^^p + C% ^ 

- E 7^ { {^-Gg - |u, - uP + - c^^) A-V . [(u, - u) V . (u, - u)]} , 

\ (a^5,) - AnGgS = -\ [aV • (J.u,)]' + ^V • (u, • Vu,) + C% l-uf^ + C^*)"^ 



^ {A [u ■ (u, - u)] - V • [(u, - u) ■ Vu + u • V (u, - u)]} 



a 

^^^^|u,; - Up - 87rGj2 Qj^ - up + SA-^V • [(u, - u) V • (u, - u)]} . (71) 



In we have shown that, to the linear order, (u^ — u) simply decays 



(66) 



(67) 



(68) 



(69) 



(70) 



Ui-uoc-, (72) 
a 



in an expanding phase. Thus, ignoring quadratic combination of (iii — u) terms, we have 

j+iv-u = --V-((5u), (73) 
a a 

iv-(u + i/u)+4^G£),5 = -^V-(u.Vu)-CW"'3 Q^,,^ + , (74) 

^ (a's) ' - iirGgS = [«V • iSu)\ + ^ V • (u ■ Vu) + C^'^^^ (^u^^p + C^^l) , (75) 
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and 



1 1 ,„ . 1 r 



^-V-u, = --V-((5»uO + - 2(^V-(u, -u)-(u, -u)-V(^ + 2«~u")l''C72 , (76) 



■ (u, + Hu,) + iirGgS = ~^\/ • (u, • Vu,) - Cf*)"'^ + (jit)\ 

1 [aH^ ' ~ ^itGqS ^~\[a^- {5,xi,)] + • (u, • Vu,) + f-<'^ + C^*)"'^ 



(77) 



+ \ {A [u • (u, - u)] - V • [(u, - u) • Vu + u • V (u, - u)]} . (78) 

Equations (|73p - (|75p coincide with the density and velocity perturbation equations of a single component medium 
0; thus, except for the contribution from gravitational waves, these equations coincide with ones in the Newtonian 
context. In the Newtonian context, Eqs. ((73|l - ([75|l without the gravitational waves were presented in [ll[ which are 
valid to fully nonlinear order. To the linear order, Eq. ([75]) was derived by Lifshitz @ in the synchronous gauge, and 
by Nariai [l9| in the comoving gauge. In the zero-pressure medium the free-falling object is also comoving, thus we 
can impose both the synchronous gauge and the comoving gauge simultaneously. In Q we compared subtle differences 
of the second-order perturbation equations in the synchronous gauge with the ones in the comoving gauge. 
If we further ignore (u^ — u) terms appearing in the pure second-order combinations, Eqs. (I76l) - (|78l) become 

+ -V-U, = -iv-(5,u,), (79) 



a a 



-V ■ (u, + Hu,) + AnGgS = -^V ■ {u, ■ Vu,) - C^*)"'' (-u^^p + cg) , 



(80) 



A^GgS = [aV • {6,xi,)] + • (u, • Vu,) -f cf^ {-uf^ + C^"'') , (81) 



which coincide with the Newtonian equations except for the contributions from the gravitational waves. In this context, 
except for the contribution from gravitational waves, the above equations coincide with ones in the Newtonian context 
even in the multi-component case. Therefore, we have shown the relativistic/Newtonian correspondence, except for 
the contributions from the gravitational waves, to the second-order perturbations in the case of multi-component, 
zero-pressure, irrotational fluids assuming a flat background. 



IV. THIRD-ORDER EQUATIONS 



We consider irrotational fluids, thus ignore all vector-type perturbations. As all three-types of perturbations are 
generally coupled in nonlinear perturbations, apparently this is an important assumption we make in this section. In 
an expanding phase, however, the linear order rotational perturbations have only decaying mode due to the angular- 
momentum conservation. Effects of rotational perturbations to the second-order perturbations are considered in our 
accompanying work in 

We take the temporal comoving gauge 

v = 0. (82) 

Together with the irrotational condition wi"-* = we have — 0. Equation ([M]) shows that this leads to = for 
general fluids. As the spatial gauge condition we take 

7 = 0. (83) 

In [l| we have shown that these gauge conditions fix the space-time gauge transformation properties completely to all 
orders in perturbations. Thus, each perturbation variable under these gauge conditions has a corresponding unique 
gauge-invariant combination, and can be equivalently regarded as gauge-invariant one to all orders in perturbations. 
We consider zero-pressure fluids, thus 

P{^) = 0, (5p(i) = EE n(i)„^. (84) 

The collective fluid quantities are nontrivial and are presented in Eq. pi)) . 
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In our previous study of second-order perturbations, summarized in a previous section, we showed that 0(1^(^)0,1^) 
terms simply correspond to pure decaying mode in expanding phase. We showed that by ignoring these terms we 
recover complete relativistic/Newtonian correspondence even in the multi-component case. Based on this observation, 
in our third order calculation in the following we will ignore Odw^j-j^p) terms. If we ignore 0{\v(i-fa\'^) terms, the 
collective fluid quantities in Eq. (|44)) give 

(85) 



and the ADM fluid quantities in Eqs. (H51) . (|46p become 

E = fi + Sfl, Ja = 0, Sal3 =0, S =0, Sap = 0, 



N ^ ait), 



Equation (|T2|) gives 



thus, 



-2i2X'"x,a(i-2^)H-lr"x'''c2. 



Using Eq. HMD, Eqs. ((n|) , ((13]) , (Ull) , (Ull) give 



Q/3l 



3iJ + iH^ + AttGq - Ac^ 



k + 2Hk - AttGqS 



,n,aX^ (1 - 2^) + + ^^•"X'''C3 



+ (^x-"' - la^'^^'^.x + [ (^x,«|, + C%) (1 - 4^) 



^ + 3i/ ) (1 + %))+%) (1 + %))^;^,) (1-2(^ + 4(^2) 



(1 - 2^) + «:<5(,) - -j,av'(,, (1 - 4^) + (1 - 4^) + '^^,5'^)X'C% 



[a (1 + %)) = (1 + (^(.,)) KW(,)a 



where we recovered c. Equation (|92p can be written as 



(1 - 2^) - 2z;['^)X'"C^S' 



(86) 
(87) 

(88) 
(89) 



(90) 



(91) 



(92) 



(93) 



Equations dMI and ^ are the same as Eqs. (20) and (21) in 3 which were derived in the single component situation. 
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Combining Eqs. (|89)) . (190)1 and Eqs. ([90)1 - ((92|) . respectively, we can derive 



1 



aH + cS^^x'" (1 - 2<^) - 2c6'^rf^cj^l\ ' - AnGgS (1 + S) 



= -k' (l + S)-^ [^oX'" (1 - 2^) + {Sk)^^ X'" - 2Af"x'''C'i*j 

-|x,.^.. - Jx- (2ci*) , - ~ < (-;Lx^V + ^^'^^0 - 4^^' 

«''^(,) + c(5(,)^„x'" (1 - 2^) - 2c<5(^)X'''C^S] ' - "^ttGqS (l + 

"^'^^ (1 + %)) + ^ [^«X'" (1 - 2^) + (<5(,)k)^^ X'" - 2«-x''^cS' 



4c 

-^X,c<P./3 



2C 



(t) 



C 



7a|/3 ck/3|7 



4c: 



(t) 



01 



A 



+ 



+2v 



"1/3 



(94) 



'S + (-4^^S + ' '^C%C^t) ] + 2i^w|. ['^ (2^ - <^(,)) + ^ {h^ - 4^)] 



-2A 



(95) 



All terms in the RHS of Eq. (|95| contain which decays to the linear order as 



V(j) = - (Uj - u) oc -, 



(96) 



see Eqs. (I60|l,(l72]l and ^5]. By setting 6(i) = 6 the LHS of Eq. ([95|) is the same as Eq. ([Ml)- By setting S^i) = 6 and 
w^^) = 0, Eq. dHSjl becomes Eq. (|M)). 

In order to complete Eq. ()95p we need equations for tp, w"^^, and C^*^ terms to the second-order perturbations. 
Equation for is in Eq. (|M)) which gives 



(97) 



to the second order. The relation between k and x can be foimd in Eq. (23) of [J. Recovering the background 



15 



curvature, from the momentum constraint equation in Eq. (I16p . we can derive 



Mt) 

■-\13 - 



to the second order. Equation for to the second order foUows from Eq. (99) in 1]. Ignoring 0{\v(^i)a\'^), thus using 



Eq. 



for a we can derive 

3(/5 = 



(99) 



In a flat background, the RHS is the second order; thus, to the hnear order we have (p — Equation for to the 
second order is presented in Eq. (43) of 



A. Pure general relativistic corrections 



Now, we assume 



= 0, 



(100) 



thus if = Q to the hnear order. Based on the apparent success in second-order perturbations, we continue to use the 
identifications made in Eq. (|63p vahd to the third order, thus 



1, 



-V • u, u = Vu, V, 



1 



{i)v 



(u,; - u) , 5 = 5v, 6i = 6(^i)y. 



(101) 



In the following we consider pure scalar-type perturbations, thus set C^*^ = 0. Contributions from the gravitational 
waves will be considered in the next section. 

We need Xv only to the second order. Equation ((98)) gives 



Xv 



(102) 



where 



X = 2(pV • u - u • V<^ -I- 2^"^^ [" ■ C^^) + "^V?] ■ 



Equations 



become 



S+-W ■u = ■ (Su) + - (2ipu - WA-'^X) ■ V6, 
a a a ' 

• (li + i7u) + 47^G^?5 = -^V • (u • Vu) - 4 [u • V iA'^X) 



-^(^u- V(V-u) + ^V 



( u • Vu — gUV • u 



S, + -V-u, = --\7 ■ {6,Ui) + - [2ipu, - V (A^^^)] • V(5, + - [2ipV ■ (u, - u) - (u, - u) • \/tp] 
a a a a 

2 1 

H — ip [6iV ■ (uj - u) + 2 (uj -u) -Vip- 2ip\/ ■ (u^ - u)] Si (u^ - u) ■ Vip, 

a a 

-V ■ (li, + Hu,) + AirGgS = ■ (u, • Vu,) - A k • V (A-^X)] +l-(u-VX + • u ) 

a \ 3 / 



-A^u.V(V.u) + ^V 



f [ u ■ Vu — 2*^^ ' ^ 



+ 2— [ifu ■ {u, - u) 



(103) 



(104) 



(105) 



(106) 



(107) 
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Equations (fTT7|) . ITTS)) coincide with Eqs. (39), (40) in Equation gives 



i[a(u,-u)]- = -iv{(l-2^) (u,-u). [u + V(A-iX)]} 



(108) 



Equation 



3 give 
aV • ((5u) - 
uVX + 
aV • ((5iUi 



a [2ipu - V (A^^^)] • V(5| - 47rGi?J = • (u ■ Vu) + ^ [u • V (A^^X)] 



3^^u.V(V.u)--iv 



(fi [ u ■ Vu — • u 



(109) 



[2ipui 



■ VS^} ~ AirGgS - • (u, • Vu,) - ^ [u, • V (A-^X) 



u-yx + -XV ■ u 



2 
3a 



2¥?U- V(V-u) 



A 



rV 



(fi [ u ■ Vu — gUV • u 



= - [2(pV • {ui - u) - (ui - u) • Vip] - 2^ [(^u • {ui - u)] 



{A [u • (u, - u)] + (V • u,) V • (u, - u)} + ^ (V^) • {V [u • (u, - u)] + (V • u,) (u, - u)} . (110) 



Equation (1991) becomes 



3a 



[-X + 2ipV • u - u • Vip] 



-u ■ V (V.^) + u—ip 
a a 



(111) 



A close examination of above equations reveals the foUowings. 

(i) Equations (|104p . (|105p . (|109p are the same as Eqs. (25), (26), (28) in [1] which are valid in a single component fluid; 
we note, however, that Odwjj-jQp) correction terms, which simply decay in an expanding phase, appear even in the 
second-order perturbations, see Eqs. ([66]) . ([67 ]) . (fTOl) above, or 

(ii) If we ignore the i-indices, Eqs. (fTOe]) . (fTOT]) . ITTOl) are identical to Eqs. (fT04| . (fT05ll . (fTog]) . respectively; we note, 
however, that the presence of Odwi^jj^p) correction terms causes differences even in the second-order perturbations, 
see Eqs. ((50| - ((55l) . 

(iii) We already showed that, to the second order, even in the case of multi-component, the general relativistic 
equations are identical to the Newtonian ones, thus having relativistic/Newtonian correspondence. The presence of 
^(b(i)Qp) correction terms may cause differences, but we have shown that these corrections are simply decaying in 
the expanding phase. 

(iv) The pure third-order correction terms in the above equations all involve ip-term to the linear order in various 
forms of convolution with the second-order terms. As we took the comoving gauge v = 0, the spatial curvature 
variable is the same as a gauge-invariant combination (fy = ip — aHv/c to the linear order. 

(v) We note that, to the linear order, is one of the well known conserved quantity in the large-scale 0, [2l|. For 
K ^ 0, but considering general A, we have (py = 0, thus 



C7(x), 



(112) 



with vanishing decaying mode (in an expanding phase) to the leading order in the large-scale expansion |T8| . 

(vi) We also note that the value of ipy in the large-scale is of the same order as the gravitational potential fluctuation 
ip^ which is again of the same order as the relative temperature fluctuations ST/T of the cosmic microwave background 
radiation (CMB). In general we have Q 



ipy ^ (p^- aHv^/c = + ciA ^V • u/c^. 



(113) 



to the linear order, where we have 
identified to the linear order in [3, Ojl, [2; 



— 5$/c^ and u = — cVw^; (5$ is the Newtonian gravitational potential 
For K = = A we have 



(114) 



The temperature anisotropy of CMB, in a flat background without the cosmological constant, gives |23l. |24 
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to the linear order; this is a part of the Sachs- Wolfe effect [H]. The COBE observations of CMB give 5T/T - lO'^ 
dl, thus 

ip^r^hX 10"^ (116) 

in the large-scale limit near horizon scale. Therefore, the pure general relativistic third-order correction terms are 
independent of the presence of the horizon scale and are smaller by factor 5 x 10~^ compared with the second-order 
relativistic/Newtonian terms due to the low level anisotropics of the cosmic microwave background radiation. 



B. Contributions from tensor- type perturbations 



Here we present set of equations describing the scalar-type perturbation equations to the third order, now including 
the contributions from the tensor-type perturbations. We continue to use the Newtonian variables identified in Eq. 

Equations (EUl-dMl) become 



(117) 



n n n ^ ' n '-^f-' 



-V • (li + Hu) + AttGqS ^ -\\7 ■ (u • Vu) - C'*)"^ 



+ ^a(3 



4[u-V(A-ix)]+l 



u-VX + -XV ■ u 



3^^u.V(V.u) + ^V 



f [ u ■ Vu — • u 



\a y [ 3a ^ ' a \ "''It/ \a I" 

1 
a 



d n n \/j n ^ 



1 



2^V ■ (u, - u) - (u, - u) ■ V<^ + 2C'i*] « - u^'f 



a 

2 1 

-I — (p [(5iV • (uj - u) -I- 2 (ui - u) • Vv? - 2(p\/ ■ {u., - u)] 5i (u^ -u) -Vip 

a a 



I" 



-V • (u, + Hu,) + AnGgS = ■ (u, • Vu,) - C^*)"'' f-u^^p + C^*^ 



A 



— [u, • V (A-i^)] + _ u • VX + -XV ■ u 



-^^u.V(V.u) + lv 



f [ u ■ Vu — gUV • u 



A 
2- 



(u, - u) + (lif - w'') 



+ 4""^^ (V • u) Cj^l + 2C'(*)"^ 



+2 (7(*)"^) [-^ (V . u) + {2C% - C^l,;) + 2cSi (^.'^ + 

Equations (fTTT)) . (fTTS)) coincide with Eqs. (39), (40) in 4]. Equation gives 



^ [a (u, - u)]' = -iv |(1 - 2(f) (Ui - u) • [u + V (A^^X)] - 2u" (wf 



/3 



(118) 



(119) 



(120) 



(121) 
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Equation (IM)) gives 

{a^^ + aV • ((5u) - a [2ipu - V (A-^X)] ■ VS - 2a(5'"M'^c3} - AnGgS 



a- 
2 

2 



V • (u • Vu) + C^*)"'^ + 



4 [u ■ V (A^^^)] -\(u-VX + IxV ■ u 



yju • V (V • u) ■ 



(yj ( u • Vu ~ ■ 



(V • u) C^*' - 2C'(*)"'' 



4 2 _ . ^(,) 1 



-2 iti"!'' + CW"'^ 



it) 

a/3|7 



2Cj*^ ( -u- 



. (u, . Vu,) ~ CW"'^ + Ci*)) - A [u, . V {A-'X)] +^i^u-VX + ^XV-u 



-A^u.V(V.u) + -iv. 



(/? ( u • Vu — gUV ■ u 



+ 4w"V^ (V • u) Cj^l + 2C(*)"'5 



,q|/3 



1 



= f [2(^V • (u, - u) - (u, - u) • V^] - 2— (^u . (u, - u) + (uf - u'^j cg 



^A 

a ' ' ' ' ' ' ' ' o? 

-\ip {A [u ■ (u, - u)] + (V • u,) V • (u, - u)} + 4 (V^) • {V [u • (u, - u)] + (V • u,) (u, - u)} 

|[u . (u. - u)]^„|^ + (V • uO - «„)|^ - aCg,^ {uj - u^) + 2a [c^}^ {u] - u^) 



Equations ([98]) , (|99l) become 

X = 2ipV • u - u ■ V^ + C^*)"'' (2uc,ip ~ aC\ 



-A-'V 



u • V (V„^) + u^A^ + u^A&^l - a^^^C^^l + 2aC(*)'3^C'g|^ + aCj^i^C'^*)'''^ 
X + 2./3V • u - u • V^ + 2C(*'"'3 [u^\p + aC^' 



3a 



A 

a 



(122) 



(123) 



^u . V (V.^) + u^-^ + u^^d:l - ^''Ctl + '2C('^''C%\, + CI^^C^'^^^ 



(124) 



(125) 



V. DISCUSSION 



In this work we have successfuUy derived pure general relativistic correction terms appearing in the third-order 
perturbations of the zero-pressure irrotational multi-component fluids in a flat background. We have ignored 0{\u — 
Uip) correction terms which simply decay in an expanding phase. Our main results are presented in Eqs. (|10ip - (|llip . 
and in Eqs. (|117p - p25p in the presence of the gravitational waves. The equations for the collective component are 
identical to the ones in the single component case. If we further ignore 0(u — u^) correction terms, which are again 
simply decaying in an expanding phase, the forms of relativistic correction terms in the individual component are the 
same as the ones in the collective component. Our results show that, even in the multi-component situation, the pure 
relativistic correction terms are smaller by a factor 



(126) 
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compared with the relativistic/Newtonian second-order terms and are independent of the presence of horizon scale. 

By taking different temporal gauge condition (hypersurface condition) we can easily introduce apparently horizon 
dependent terms with arbitrarily huge amplitudes. The exact relativistic/Newtonian correspondence to the second- 
order perturbations was available essentially due to our proper choice of gauge conditions and correct identifications of 
relativistic variables with the Newtonian ones. In our third-order extension we have assumed the same identification 
holds even in the third-order perturbations, which might not be necessarily the unique choice. However, the properties 
(the smallness of the amplitudes and independence from the horizon scale) of our third-order correction terms assure 
that our choice of the gauge and identifications are very likely to be the correct and best ones even to the third order. 

Our results may have practically important implications in currently favored cosmological pursuits by assuring the 
use of Newtonian physics in the large-scale nonlinear processes which often involve two-component zero-pressure fluids 
(say, dust and cold dark matter). As we have shown the exact relativistic/Newtonian correspondence to the second 
order and small horizon-independent third-order correction terms, it is now more secure to use Newtonian physics 
near (and even beyond) the horizon scale which is indeed a noticeable trend in current cosmological simulations [26| . 
At this point it might be also worth emphasizing the effects of rotational perturbations. In Q we have shown that 
the rotational perturbations to the second-order also have relativistic/Newtonian correspondence in the small-scale 
(sub-horizon scale) limit. As the numerical simulations naturally include the rotational mode, this might be another 
good news to the cosmology community based on Newtonian physics. 

Although we have estimated that the third-order pure general relativistic correction terms are small in the current 
large-scale structures, it would be still interesting to see whether cosmology could reach a stage where ~ 5 x 10^^- 
factor smaller correction terms could have noticeable effect on the large-scale structure formation processes. In either 
case the third-order terms we have presented in this work are the first non-vanishing general relativistic correction 
terms in Newtonian nonlinear equations, and may have historical as well as practical importance in cosmology. 
Analytic studies in the Newtonian context of a single component zero-pressure, irrotational fluid have been actively 
investigated in [27j . Whether our general relativistic corrections appearing in the third order (both in the single- 
component and multi-component cases), and corrections appearing even in the second order by effects of pressure, 
rotation, and background curvature, will have noticeable roles in the large-scale structure formation process are left 
for future investigations. 
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